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High-amplitude acoustic waves in thermoacoustic devices are responsible for the generation of a second¬ 
ary flow called acoustic streaming. Superimposed on the oscillating flow, this secondary flow is an impor¬ 
tant source of energy dissipation. To remove acoustic streaming would result in substantial 
improvements in the energy performance of thermoacoustic devices. Understanding the control param¬ 
eters and mechanisms of streaming generation is essential for controlling acoustic streaming. 

In this paper, streaming sources are investigated in a pressurized thermoacoustic prime mover. The 
device was designed and built at the LIMSI (Laboratoire d'lnformatique pour la Mecanique et les Sciences 
de l’lngenieur) to investigate acoustic streaming. To calculate the theoretical velocity in the thermoacou¬ 
stic device, the usual modelings were extended to take into account the cross section variation. Measure¬ 
ments of acoustic streaming velocity performed by Laser Doppler Velocimetry are compared to the 
theoretical results. Finally, the good agreements obtained led to the validation of the theory. The vali¬ 
dated theory was then used to study the mechanisms of streaming generation in order to determine their 
influences in the thermoacoustic prime mover and therefore reducing the secondary flow. Thereafter, the 
investigation was extended to four channels with different widths and according to results of the litera¬ 
ture, the proportion of outer streaming compared to inner streaming decreases when the channel width 
decrease. The investigation showed that the mechanisms of acoustic streaming generation differ for large 
and narrow channels. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

In the energy transition century, thermoacoustic systems have 
been receiving a renewed interest. Thermoacoustic devices are a 
part of the cleanest energy conversion systems for the environ¬ 
ment because they use non-polluting gases [1,2]. As thermoacou¬ 
stic systems have no mechanical moving component, they are 
also highly reliable and have a long lifetime. To keep these reliable 
and environmental properties, pulse tube refrigerators could be 
driven by thermoacoustic engines [3-5], Despite its strengths, en¬ 
ergy efficiency of thermoacoustic devices is still low. Therefore, 
increasing the performance of these systems is an important issue 
in the thermoacoustic community [6-8]. The high mean pressure 
amplitude in these machines, necessary for their functioning, is 
responsible for the appearance of a steady mass flow of the second 
order. Usually called streaming, it is superimposed on the first 
order oscillating acoustic mass flow. These dissipating energy 
phenomena have their origin in the nonlinear nature of the propa¬ 
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gation. From energy considerations and despite their low level, 
these second order phenomena involve heat transfer to the wall, 
which is an undesirable loss mechanism. As acoustic streaming 
has been identified as an important source of energy dissipation 
in thermoacoustic devices [9-11], a better comprehension of this 
phenomenon is necessary for improving their performance. 

By introducing the Reynolds stress tensor for acoustic waves, 
Lighthill was the first to attempt to describe the mechanisms at 
the origin of acoustic streaming [12]. The Reynolds tensor, well 
known in the theory of turbulence, has been extended to the the¬ 
ory of sound and can be written in the form, pUiUj, where u, are 
the fluctuating velocities in the sound and ■ indicates the time 
average of the quantity underneath. Considering an area S oriented 
in the x* direction, pitjUjS represents the time averaged transport of 
momentum in the x, direction across S. As the Reynolds tensor acts 
as a stress, the difference between the force acting on the two 
opposite sides of a small block of fluid separated by a distance h 
generates a net force -hS[dpu t Uj /9x f [. Per unit of volume, the ith 
component of the force acting on the fluid is 

F I = -"P. (1) 





















R. Paridaens et al. / Cryogenics 58 (2013) 78-84 


79 


Nomenclature 



Co 

speed of sound (m s~') 

Ua\ 

first order of the acoustic axial velocity (m s _1 ) 

drh 

mass accumulation per unit of volume (kg nr 3 s _1 ) 

Ui 

ith component of the particle’s velocity (m s _1 ) 

diW 

mass accumulation in an elementary volume (kg s _1 ) 

UL 

average mass transport velocity (m s -1 ) 

dr 

elementary volume (m 3 ) 


streaming axial velocity (m s _1 ) 

Fr 

volumic force induced by the spatial variation of the 

Va 

acoustic transversal velocity (m s~') 


Reynolds stress tensor (N nrr 3 ) 

Va\ 

first order of the acoustic transversal velocity (m s _1 ) 

Fr, 

ith component of the volumic force due to the spatial 

V 

particle’s velocity (m s _1 ) 


variation of the Reynolds stress tensor (N nr 3 ) 

v a 

acoustic velocity (m s~') 

Fru 

first term of the volumic force due to the spatial varia- 

V„1 

first order of the acoustic velocity (m s~') 


tion of the Reynolds stress tensor (N nr 3 ) 

V„2 

second order of the acoustic velocity (m s _1 ) 

Frv 

second term of the volumic force due to the spatial var- 

v m 

streaming velocity (m s~') 


iation of the Reynolds stress tensor (N nr 3 ) 


ith coordinate (m) 

Fr 

volumic force induced by the spatial variation of the vis- 

X 

axial coordinate (m) 


cous stress tensor (N nr 3 ) 

y 

transversal coordinate (m) 

Fpi 

ith component of the volumic force due to the spatial 

a 

coefficient depending on the fluid 


variation of the viscous stress tensor (N nr 3 ) 

<>ik 

Kronecker delta 

k 

wave number (m 1 ) 

S v 

thickness of the boundary layer (m) 

P 

total pressure of the fluid (Pa) 

n 

dimensionless transversal coordinate 

Pa 

acoustic pressure (Pa) 

y 

ratio of heat capacity at constant pressure and constant 

Pat 

first order of the acoustic pressure (Pa) 


volume 

Pal 

second order of the acoustic pressure (Pa) 

p 

shear viscosity (Pa s) 

Pm 

streaming pressure (Pa) 

Pat 

first order of the acoustic shear viscosity (Pa s) 

Po 

mean pressure (Pa) 

Pb 

bulk viscosity (Pa s) 

R 

width of the channel (m) 

Po 

mean shear viscosity (Pa s) 

Re[] 

real part 

CO 

angular frequency (s _1 ) 

s 

transverse section (m 2 ) 

% 

complex functions dependent on the thermophysical 

Sa 

mass flow rate due to the propagation of the acoustic 


properties of the fluid and geometric parameters, 


wave (kg nr 3 s _1 ) 


i = 1,2,3 

Si 

elementary surface of the volume oriented toward x,- 

<l>i 

contributions of the streaming sources, 


(m 2 ) 


i = 1,2,3,4 (m s _1 ) 

T 

total temperature of the fluid (K) 

p 

total density of the fluid (kg nr 3 ) 

T a 

acoustic temperature (K) 

Pa 

acoustic density (kg nr 3 ) 

Tat 

first order of the acoustic temperature (I<) 

Pai 

first order of the acoustic density (kg nr 3 ) 

To 

mean temperature (K) 

Pal 

second order of the acoustic density (kg nr 3 ) 

T, 

relative contributions of the streaming sources (%) 

Pm 

streaming density (kg nr 3 ) 


i = 1,2,3,4 

Po 

mean density (kg nr 3 ) 

U 

magnitude of the acoustic velocity (m s~ ] ) 

G 

Prandtl number 

U a 

acoustic axial velocity (m s~') 

<?ik 

viscous stress tensor (Pa) 

Uai 

ith component of the acoustic axial velocity (m s~') 




This force contributes to acoustic streaming generation [13], 
Another well known mechanism at the origin of the secondary 
flow is the effect of viscosity. Considering an oscillating parcel of 
gas located at about the penetration depth from the wall, Olson 
and Swift described this phenomenon [14], Due to thermal contact 
with the wall and the phasing between oscillatory pressure and 
motion, the gas between the parcel and the wall has a varying tem¬ 
perature during the oscillation. The parcel experiences differing 
amounts of viscous drag during its motion because the viscosity 
depends on the temperature. After a full cycle, the parcel does 
not return to its initial point. Bailliet et al. have developed a model 
describing streaming for standing waves and an arbitrary channel 
width, in which the temperature dependence of the viscosity is in¬ 
cluded [15]. In that model, the volumic force due to viscosity on 
the axial component x is 

F /( x = div, ; (^grad(u)), (2) 

where x is the axial component and i( is the transversal one. u is the 
x-component of the acoustic velocity and p is the shear viscosity. 
The time average of this force is non-zero because the viscosity is 
time-dependent. The volumic force generated by the viscosity can 
be written in the form 


As for the Reynolds stress, the volumic force at the origin of 
streaming is due to a spatial variation of a stress tensor. The 
acoustic mass flow can also be a cause of streaming generation. 
The Reynolds stress, the viscosity, and the acoustic mass flow, 
identified as the phenomena at the origin of acoustic streaming 
will be quantified in a pressurised thermoacoustic primemover. 
Finally, the investigation on streaming sources will be extended 
to four channels with different widths. 

2. Mechanisms of acoustic streaming generation 

This section focuses on the mechanisms of acoustic streaming 
generation. The streaming equations are presented, discussed 
and solved. To investigate nonlinear effects in thermoacoustic de¬ 
vices, some assumptions are made. A secondary flow is assumed 
to be superimposed on the oscillating flow. The time and space 
dependent variables can be written in such a way as to take into 
account the variables in the absence of oscillations {p 0 .p 0 ), acous¬ 
tic variables (p a ,p a ,v a ), and secondary flow variables [p m ,p m , v m ): 
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P = P 0 + Pa + Pm> P = Po + Pa+Pnv V=V a +V m . (4) 

Each acoustic variable could be developed asymptotically: 

Pa = Pal+Pa2 +Pa = Pal +Pa2 + ■ ■ ■ , V a = V al +\ a2 + .... (5) 


The acoustic Mach number is assumed to be much smaller than 
1. The streaming is assumed to be stationary and the Mach number 
associated to the secondary flow is less than 0.3. Therefore, the 
secondary flow is expected to be incompressible. The Reynolds 
number of the streaming is significantly less than 1 , so the advec- 
tion term is negligible compared to the viscosity term. With this 
assumption, the corresponding secondary flow is called slow 
streaming. In this paper, the investigations focus on the Rayleigh 
streaming, and the mass flow through any cross section is zero in 
the steady-state operation. Under these assumptions, the equa¬ 
tions governing the streaming at second order are [15]: 


P o 


du m 
dx + 


dv m 

dy. 


= S a , 


( 6 ) 


0 = - 


9Pn, 
dx 


+ Po 


9%r 

dy 2 


+ Fr + Pi 


n, 


( 7 ) 


^~ 0 . 
dy ’ 

where 


Sa=- 


9 Pa 1 Uni 
dx 


9P a i Val 
dy 


( 8 ) 

( 9 ) 


Fr 


' dp 0 u aj u a j dp 0 v a iu a , 
dx dy 


( 10 ) 


h 


d_ 

9y 


PaV 


9Ua\ 

dy 


( 11 ) 


Eqs. ( 6 )-( 8 ) are the mass and momentum conservation equa¬ 
tions governing the streaming, u m and v m are, respectively, the ax¬ 
ial and the transversal velocity of the acoustic streaming, p m is the 
additional pressure generated by the streaming, u ai and v a i are, 
respectively, the axial and the transversal acoustic velocity. As 
the shear viscosity depends on the temperature [16], it can be writ¬ 
ten as 

p oc T p . (12) 

P is a coefficient depending on the fluid and /( al can be expressed as 
a function of T„i 

Pal = PoP^f -1 ( 13 ) 

where T a i is the additional temperature due to the acoustic wave. 
S a ,F R , and are the acoustic streaming source terms, which are 
only governed by the acoustic field. F R and F fl are volumic forces 
due to, respectively, the spatial variation of the Reynolds tensor 
and that of the viscous stress tensor. F R is composed of two terms: 


9p 0 UalUai 

dx ’ 

(14) 

dp Q V a \U a \ 

dy 

(15) 


The term S a is a mass flow rate induced by the propagation of 
the acoustic wave. To understand the meaning of this term, it is 


if > (\ !l UlPl)[z L ) 


> 


(pQU ai Si)( Xi 

+dx.i ) 


> 


x i Xi + dxi 

Fig. 1. Mass balance in an elementary volume dr = dx, dx, dx„. 


assumed that there is no secondary flow. With this assumption, 
the density and the velocity are 

P = Po + Pai u = u al . (16) 


The mass flow accumulated per unit of time inside an elemen¬ 
tary volume dr (Fig. 1) is 

dJW = [p^Uai(Xi) - J)^Uai(Xi + ( 17) 

S, = dXjdXfc |. k/ . is an elementary surface of the volume dr. Per unit 
of volume, this mass accumulation becomes 


_dM_ dp aX u ai _ 

di dXi 


(18) 


Without average flow, S a represents a mass accumulation per unit 
of volume. As there is no accumulation, this mass is evacuated by 
generating a secondary flow, referred to as acoustic streaming. 

To calculate the axial streaming velocity in a thermoacoustic 
device, the modeling of Bailiiet et al. [15] was extended by incorpo¬ 
rating the section variation. The solution of Eqs. ( 6 )—(11) to second 
order gives an expression for the axial streaming velocity: 

Um = l^l + 4*2 + <(>3 + <^ 4 , (19) 


with 


3 /'• 

h = 4 ^ o (' 7 2 ~ 1 ) J MZbl, 


( 20 ) 


<l>2= — 
Po 

_3 

4 


<t> 3 = — 
V 0 


n 

(i ->f) /' f [ 

J -1 J 1 Jo 


9(P 0 UalU a i) 

dx 

1 nY\ rt]' 


difdrj' 

d(p 0 UaiU a i ; 

dx 


dtfdifdi] 


( 21 ) 


, P 

<P4 = ~r~ 

J o 


rl q r< rn 

J U al v aX dif - -(1 - if) J ^ j u 0l Va^difdi] 


( 22 ) 


(23) 


The contribution of the mechanisms mentioned previously are 
represented by </> ( in Eq. (19): </;,, r/> 2 , </) 3 , and <j> 4 are, respectively, 
the contributions of the term S, F Ru ,F Rl „ and F ;i to the axial streaming 
velocity, and i , i = y/R is the normalized transversal dimension. The 
centerline and the wall are respectively located at rj = 0 and rj = 1 . 

The acoustic variables are expressed in function of the acoustic 
pressure in a channel of non-constant section. The conservation 
equations of mass, momentum and energy and the state equation, 
solved in 2D Cartesian coordinates, yield the first order velocities, 
temperature and density: 
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Wfll 


1 p dp a i 

icop 0 11 dx 


f _?~1t r „ 1 f>/t oFtf dT o 9p al 

1 “1 — „ ,.2 ‘OrytPai , 1)2 , 


Po c o 


co 2 p 0 1 - a dx dx 


y - (7 - 1 )F,,t „ , 1 fy - crF, dT 0 dp aA 

Pal —-3-Pal + 


icoR 
Va\ =TTT2 
Po c o 


r 0 ct ) 2 1 - a dx dx 

7 + (1 - 7 ) F t 


(24) 

(25) 

(26) 


yy + (i - 7)( f / - <W + (</>- '7) J 


Pal 


iRb 2 F) p(F - \) + t y dR dp a , 


fflPo V 

- — (, 
mp 0 V 


R <9x 9x 


b 2 (l - F) ~ J) + ^ + r]F, 


/? + 1 9F 0 9p a i 
2 T„ ~dx~dx 


+ 


i'R 




where 

b= (l- 0 K 


<5v = 



, sinh(bw) , sinh (\fobif) 

<P = L mr L,.,' , <Pt = ] --- 


bcosh(b) ’ 


byircosh (bv/ff) 


(27) 

(28) 

(29) 


r cosh(bp) . tanh (i/ffb) tanh(b) 

= Ft=1 -’ F=1 -b~ 

(30) 


cr is the Prandtl number and y is the ratio of heat capacity at con¬ 
stant pressure and constant volume. The expression of the axial 
velocity, the temperature and the density (24)-(26) are not modi¬ 
fied with the consideration of a section variation [10,15,17], 

The Eqs. (19)-(30) lead to an expression for the streaming axial 
velocity in terms of the control parameters: 


u m = 


Po Po c o 


Re 


R 


p 0 co z p 0 


'Pi Pal 


RepFs] 


(9PaA 

*• 

+ 

V dx) 


1 dTo 

9Pa 1 

T 0 dx 

dx 


__RePP l! — 

/( 0 cu 2 p 0 K [ 2 'Rdx 


dPa i 


dx 


(31) 


where V I / 1 , V P 2 and ¥3 are complex functions depending on the ther¬ 
mo-physical properties of the fluid (its viscosity, density ...) the fre¬ 
quency and geometric parameters. Eq. (31) differs from the 
modeling of Bailliet et al. by the second term [15]. Indeed, this term 
takes into account the section variation. Eq. (31) associated with the 
knowledge of acoustic pressure allows to calculate the axial stream¬ 
ing velocity. The acoustic pressure could be calculated numerically 
by the software DeltaEC [18], so that no assumption is made on the 
type of wave, standing or progressive. Therefore, the combination of 
a numerical and an analytical modeling allows to determinate the 
axial velocity of slow streaming in every type of thermoacoustic 
device. Eq. (31) also allows highlighting three parameters which 
control the acoustic streaming: the spatial variation of the magni¬ 
tude of the wave, the section variation, and the temperature 
gradient, As the spatial variation of the magnitude of the wave is 
the major control parameter, no acoustic streaming is generated if 
the magnitude wave does not spatially change. For example, 
without viscous dissipation, the propagation of progressive waves 
does not generate acoustic streaming. As investigated previously 
by Olson and Swift, acoustic streaming could be locally suppressed 
by setting an adequate section variation. The theoretical values of 
acoustic streaming are obtained by Eq. (31). 


3. Verification test 


Before studying the mechanism of streaming generation, the 
theoretical fields are compared to the experimental values in a 
thermoacoustic prime mover for acoustic pressure, axial acoustic 
velocity and axial streaming velocity. A thermoacoustic-Stirling 
prime mover illustrated in Fig. 2 has been designed and built espe¬ 
cially to investigate acoustic streaming. The experimental appara¬ 
tus is composed of a closed loop and a straight resonance tube. This 
energy conversion system does not have any moving part, which 
increases its reliability [19]. Electric supply is used to bring heat 
to the hot heat exchanger and water stream passes through the 
cold heat exchanger. The 4.25-m-long resonator associated with 
the closed loop allows the system to work at the frequency of 
22 Hz. The resonator is composed of three parts: a 56-mm-diame- 
ter tube connected by a tapered tube to a larger tube of 162 mm 
diameter. The system using nitrogen can operate up to a mean 
pressure of 30 bars. Kistler piezoelectric pressure sensors were 
flush mounted along the traveling wave engine allowing the 
acoustic pressure measurements. They are located from Cpl to 
Cp 8 positions. The measurements of acoustic and streaming veloc¬ 
ity are performed by Laser Doppler Velocimetry (LDV). The system 
operates with a 600-nm-length-wave laser which delivers a power 
of 35 mW. The focus of the laser beams into the pressurised 
thermoacoustic prime mover is made through an optical device 
designed and built to resist to a pressure up to 40 bars. Two optical 
devices are placed on the thermoacoustic prime mover at the 
positions x = 2.25 m and x = 3.52 m. The acoustic flow is visual¬ 
ised by injecting particles of magnesium oxide into the system. 

After setting the system at the desired pressure and supplying 
the heat exchanger, it is necessary to wait for the wave to reach 
its steady state. The time evolution of acoustic pressure repre¬ 
sented on Fig. 3 shows that the transient state lasts approximately 
90 min. t = 0 represents the moment of setting the supplier. The 
steady state is considered reached when the acoustic pressure 
has a variation less than 2 % on a period of 10 min. 

The acoustic pressure distribution along the resonator is repre¬ 
sented on Fig. 4. Experimental data are compared to the theoret¬ 
ical results for 3 different mean pressures (a) and 3 different heat 
powers (b). The theoretical acoustic pressure is numerically calcu¬ 
lated by DeltaEC. Although the channel has a circular section, the 
calculations are performed by a 2D-Cartesian model. Therefore, 
the surface curvature is assumed not to affect the studied phe¬ 
nomena. This assumption has already been used in the literature 
[20,21], The theoretical results match the experimental ones with 
a relative error less than 6%, which validates the numerical 
modeling. The pressure nodes, depending only on the geometric 
configuration, are located at the position x = 2.23 m regardless 
the heat powers and the mean pressures. As the acoustic and 
streaming variables are expressed in function of the acoustic 
pressure and its derivative, the determination of the theoretical 
acoustic pressure field allows to calculate the acoustic and 
streaming velocities. 

The theoretical velocities are represented by the lines in Fig. 5 
for a mean pressure of 10 bars and three different heat powers: 
165 W, 190 W and 210 W. They are compared to the experimental 
data performed by LDV and good agreements are observed. 
Concerning the acoustic streaming, as mentioned in the literatures 
[12], two different velocities could be used. The first velocity, 
expressed in Eulerian coordinates, is u m obtained from Eq. (19). 
The second one, expressed in Lagrangian coordinates, is called 
the average mass transport velocity u L . The relation between the 
Eulerian and the Lagrangian velocity is 


u L = u m + 


Pa 1 U a1 
Po 


(32) 
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Thermal buffer 
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heat exchanger 




i 
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Fig. 2. Experimental setup. 


x 10 5 



Time (min) 


Fig. 3. Time evolution of the acoustic pressure. 


As the streaming velocity obtained with LDV techniques have to 
be compared with the Eulerian streaming velocity, u m from Eq. (31) 
is plotted versus ;; in Fig. 5 for x = 3.52 m. The experimental data 
performed by LDV are represented in dot points. A good agreement 
is found between experimental and calculated results. This config¬ 
uration allows to compare the order of the acoustic velocity with 
the streaming velocity; it is observed that the streaming velocity 
is two orders smaller than acoustic velocity. Similar results, not 
represented in this paper, are also obtained for all combinations 
of a mean pressure of 10, 14 and 18 bars with a heat power of 
165 W, 190 W et 210 W. 


4. Discussion 

As the acoustic and streaming fields have been obtained and 
validated by experimental data, the influence of acoustic streaming 
sources can now be investigated. Thus, the different contributions 
are plotted versus 17 for x = 3.52 m on Fig. 6 for a mean pressure of 




Fig. 4. Acoustic pressure distributions for a heat power of 210 W and three mean 
pressures 10, 14,18 bars (a) and for a mean pressure of 10 bars and 3 heat powers 
165 W, 190 W, 210 W (b). (—) Theory, (■) experimental values. 


10 bars and a heat power of 210 W. The absolute contributions (a) 
are calculated using Eqs. (20)-(23) and the relative ones (b) are ob¬ 
tained by the following equation: 

t. =_l^il_ (331 

' M + lfcl + lfcl + W 1 1 

Outside the boundary layer, the relative contributions maintain the 
same values whatever the location. In this configuration with a con¬ 
tribution equal to 92%, the Reynolds stress are the main mechanism 
of the streaming velocity generation. This result is in agreement 
with the theory of Rayleigh. Because Rayleigh neglected the tem¬ 
perature dependence of viscosity and considered that the streaming 
is divergence free, the streaming calculated by Rayleigh is only due 
to Reynolds stress <j> 2 and (j> 3 : 

tfmRayleigh ~ 02 T 03 (34) 
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Fig. 7. Absolute (a) and relative (b) contributions of streaming velocity versus x for 

'1 = 0. ( 1 ) 0,, ( ) 02, ( ) 0 3 , (--) 04- 





Fig. 5. Acoustic and streaming velocity profiles for a mean pressure of 10 bars and 
three differents heat powers: 165 W, 190 W and 210 W. (—) Theory, (■) LDV 
measurements. 


As his modeling predicts in a meaningful way the streaming 
velocity in a wide channel, it is coherent to find that the major 
source of streaming generation is due to Reynolds stress in our 
configuration. In our case, for R/S v =360,Ti,T 2 ,r 3 , and I 4 make, 
outside of the boundary layer, a contribution of, respectively, 1%, 
55%, 37%, and 7%. Similar values are also obtained for mean pres¬ 
sures of 10, 14 and 18 bars and heat powers of 165 W, 190 W 
and 210 W. The values of the relative contributions seem not 
depending on the magnitude of the acoustic wave. 

Also investigated versus x, the absolute (a) and relative (b) con¬ 
tributions are plotted in Fig. 7 for i] = 0, a mean pressure of 10 bars 
and a heat power of 210 W. The relative contributions maintain the 
same values whatever the location. The velocities vanish at the po¬ 
sition x = 0, corresponding to the condition without penetration 
into the wall, and at the position x = 2.23 m, corresponding to 
the acoustic pressure node. The axial acoustic velocity is oriented 
towards the increasing x from 0 to 2.23 m to change direction be¬ 
tween 2.23 m and 4.41 m. Therefore, in the whole resonator, the 
streaming velocity for ;; = 0 is always oriented towards the acous¬ 
tic pressure node. This result is in agreement with the investigation 
of Hamilton et al. [22]. As the acoustic pressure has a higher ampli¬ 
tude from 2.23 m to 4.41 m than that from 0 to 2.23 m (Fig. 4), sim¬ 
ilar result is observed for streaming velocity. The influence of the 



<|) (m/s) 



100 


Fig. 6. Contributions of streaming velocity versus 77 for x = 1.62 m (a), and versus x 
for ii = 0 (b). ( ) 0„ (—) 0 2 , ( ) 03, (-) 0 4 . 






Fig. 8. Streaming velocity and its contributions versus rj for different channel 
widths. (-) U m , (,■,) 0,, (—) 02, ( ) 03, (-) 0 4 . 
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section variation on the streaming velocity, represented by the sec¬ 
ond term of Eq. (31 ), could not be seen on Fig. 7 as in this case the 
tapered tube is located at the vicinity of the acoustic pressure node 
x = 2.23 m. The Rayleigh vortex is located at this position and the 
axial streaming velocity is almost zero. 

In order to understand the physical interpretation of the differ¬ 
ent contributions, the present study was carried further: the influ¬ 
ence of the channel width on the acoustic streaming was analysed. 
The study was performed in a 1 m long rectangular resonator sub¬ 
jected to a standing wave. In the resonator, the mean pressure is 
10 s Pa and the mean temperature is 293 K. The wave frequency 
is 170 Hz, a resonant frequency of the device functioning with 
air. The pressure magnitude is chosen to be 10 5 Pa. The investiga¬ 
tions were made versus the ratio R/5 V as this is a pertinent param¬ 
eter with which to quantify the channel width. 

In Fig. 8, the axial acoustic streaming (u m ) and its contributions 
(0,, <t> 2 , </> 3 , </> 4 ) are plotted versus ij = y/R for R/5 V equal to 100, 20, 
10, and 1. As inner streaming is delimited by a succession of a max¬ 
imum and a minimum of the velocity, an increase is observed 
when R/S r decreases from 100 to 10. The inner streaming is due 
to <j > 2 as the sequence of a maximum then a minimum appears only 
for the first component of the Reynolds stress. Unlike inner stream¬ 
ing, outer streaming is characterized by a sequence of a local min¬ 
imum then a local maximum. This sequence appears for tj> 2 and <j> 3 . 
Outer streaming is due to the Reynolds stress. Because both inner 
and outer streamings are due to the Reynolds stress, they naturally 
disappear for R/5 V = 1 when becomes predominant (Fig. 8 (d)). 

5. Conclusion 

Acoustic streaming has been investigated in a pressurised ther¬ 
moacoustic prime mover. To calculate acoustic streaming, the 
equations were solved to take into account the section variation. 
Then, measurements of the streaming velocity were made and 
compared with the theoretical values. The experimental values 
fit the model which allowed the validation of the theory. Then, 
the mechanisms of streaming generation were investigated. Three 
phenomena causing the acoustic streaming have been identified: 
acoustic mass flow, Reynolds stress, and viscosity. Their physical 
interpretations have been given and their contributions have been 
evaluated. The mechanism of mass flow is found to be the main 
contribution of acoustic streaming velocity in narrow ducts 
whereas the major contribution in large channel is due to Reynolds 
stress. When the channel width is decreasing, the proportion of in¬ 
ner streaming is increasing. Because inner and outer streaming is 
due to Reynolds stress, at a certain limit when the channel is too 
narrow, there is no inner and outer streaming. This investigation 
and the understanding of the origins of streamings are an impor¬ 
tant step in improving the energetic performance of thermoacou¬ 
stic systems. 
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